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Abstract. Most of the inspiralling compact binaries are expected to be 
circularized by the time their gravitational-wave signals enter the frequency band 
of ground-based detectors such as LIGO or VIRGO. However, it is not excluded 
that some of these binaries might still possess a significant eccentricity at a few 
tens of hertz. Despite this possibility, current search pipelines - based on matched 
filtering techniques - consider only non-eccentric templates. The effect of such 
an approximation on the loss of signal-to-noise ratio (SNR) has been investigated 
by Martel and Poisson, Phys. Rev. D. 60 124008, in the context of initial LIGO 
detector. They ascertained that non-eccentric templates will be successful at 
detecting eccentric signals. We revisit their work by incorporating current and 
future ground-based detectors and precisely quantify the exact loss of SNR. 

In order to be more faithful to an actual search, we maximized the SNR over 
a template bank, whose minimal match is set to 95%. For initial LIGO detector, 
we claim that the initial eccentricity does not need to be taken into account in 
our searches for any system with total mass M G [2 — 45] if eo 0.05 because 
the loss of SNR (about 5%) is consistent with the discreteness of the template 
bank. Similarly this statement is also true for systems with M G [6 — 35] and 
eo ^ 0.10. However, by neglecting the eccentricity in our searches, significant loss 
of detection (larger than 10%) may arise as soon as eo > 0.05 for neutron-star 
binaries. We also provide exhaustive results for VIRGO, advanced LIGO and 
Einstein Telescope detectors. It is worth noticing that for Einstein Telescope, 
neutron star binaries with eo > 0.02 lead to a 10% loss of detection. 



PACS numbers: 02.70.-c, 07.05.Kf, 95.85.Sz, 97.80.-d 



1. Introduction 

Inspiralling compact binaries are one of the most promising source of gravitational- 
wave signals for interferometric ground-based detectors such as LIGO [1] and 
VIRGO [2] . The phase and amplitude of the signal strongly depend on the masses of 
the two compact objects, the spin effects (neglected in this paper) and eccentricity. 
However, it is known that most of the inspiralling compact binaries would have 
circularized by the time they enter the lower cut-off frequency of the detectors {Fl 
hereafter, which is about 20 Hz to 40Hz). Indeed, the eccentricity of an isolated 
compact binary evolving under the effect of gravitational radiation reaction is reduced 
by a factor of three when the semi-major axis is halved [3]. For instance, even though 
the orbital eccentricity of the Hulse and Taylor binary pulsar is eo = 0.617 (from the 
last measurements published in [4]), it would decrease to 10~^ by the time the orbital 
period of the binary reaches 0.02 s [5]. Therefore, any isolated binary star formed 
by stellar evolution would have becomes circularized by the time its signal become 
measurable by ground-based detectors. 

Nevertheless, in addition to the usual route to binary star formation, there 
exists a different formation scenario in densely populated regions (e.g., in globular 
clusters) where binary-binary interactions can produce hierarchical triplets of black 
holes. The inner binaries of these triplets (undergoing Kozai oscillations) can merge 
under gravitational radiation reaction with large eccentricity [6]. Consequently, even 
at the lower cut-off frequency of ground-based detectors, there are compact binaries 
with significant eccentricity. As predicted in [7], about 30% of the binaries will posses 
eccentricities larger than 0.10 when their emitted gravitational- wave signals reach a 
frequency of 10 Hz. There is also semi- analytical description of the final stages of 
mergers of black-hole neutron-star(BHNS) binaries showing that they may transit to 
inspiralling compact binaries in eccentric orbits [8]. 

More recently, in [9] , the authors study the case of stellar mass black holes (BHs) 
that segregate and form a steep density cusp around supermassive black holes in 
galactic nuclei. They found that BH binaries that form this way in galactic nuclei 
have an expected rate of coalescence events detectable by Advanced LIGO about 
1 — 1000 per year, depending on the initial mass function of stars in galactic nuclei 
and the mass of the massive BHs. The authors also claimed that such BH binaries 
have significant eccentricities as they enter the LIGO band (90% with e > 0.9), and 
hence they are distinguishable from other circularized binaries. They also show that 
mergers of eccentric binaries can be detected to larger distances and greater BH masses 
than circular mergers, up to (~ 7OOM0). 

Despite the fact that binaries may have non- negligible eccentricities, current 
searches for inspiralling compact binaries neglects this parameter. For instance, the 
matched filtering technique that is deployed to analyze scientific runs of LIGO data 
considers templates with circularized orbits only [10]. Yet, matched filtering is sub- 
optimal if the template is not a faithful representation of the signal that is searched for. 
In a previous work, Martel and Poisson [11] studied the loss of SNR when signals from 
eccentric binaries are filtered with templates based on circularized binaries. They 
performed this study in the context of initial LIGO. For instance, in the case of a 
(1.4 + 1.4)M0 binary with an initial eccentricity of 0.05, they found that only a 2.4% 
loss of SNR is expected. This result combined with the idea that eccentricities are 
small when a gravitational-wave signal enters a detector's band, led to searches where 
eccentricity is systematically neglected. 
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In Ref. [11], the authors considered a continuous space of templates to compute 
the loss of SNR. In this paper, we consider a more realistic situation by using a 
discrete grid of templates, which is identical to the one used in LIGO searches (where 
eccentricity is neglected). We also extend our target sources not only to binary 
neutron star (BNS) and low mass binary black hole (BBH) but also to large mass 
BBH and BHNS. Our goal is to provide a precise description of the loss of SNR when 
waveforms from eccentric inspiralling compact binaries are searched with waveforms 
from circularized binaries. 

In Sec. 2, we describe the relevant features of the inspiraling compact binaries. 
We also emphasize the differences between our waveforms and those obtained in [11]. 
In Sec. 3.2, we describe the simulation protocol and present exhaustive results showing 
the loss of SNR for different detectors, namely initial LIGO, advanced LIGO, Einstein 
Telescope and VIRGO detectors. Finally, in Sec. 4, we discuss the implications of our 
results for ground-based detectors and future directions for this work. 



2. The eccentric model 

2.1. eccentric waveform calculation 

The evolution of the orbital elements of an eccentric binary was first calculated by 
Peters and Mathews [3, 12]. The explicit expressions of the gravitational- wave signals 
were first derived in [13], in the context of spacecraft Doppler detection of gravitational 
waves. Although eccentric waveforms have been derived to a higher Post-Newtonian 
order [14], we restrict ourselves to the case of Newtonian phase order. We use the 
same notation as in [11]. Let us start from the orbital radius r that is given by 

r= (1) 

1 + e cos (j) ' 

where p is the semi-latus rectum, M is the total mass, e is the eccentricity, and (j) is 
the orbital phase. The semi-major axis is related to p and e by the relation 

pM 
1 — 

Note that the orbital period is defined by 



(2) 



The decay of orbital elements p and e are calculated in the quadrupole 
approximation [3]. Together with the orbital phase evolution, we have three ordinary 
differential equations (ODE) that can be used to describe the evolution of the system 
and form the starting point to determine the eccentric waveforms: 

\2 



d(j) (l + ecos0)^ 



(4) 



dt p2M 

dp _ 64 n f , . 7 

dt~ 5 p3 

de _ 304 M , 121^2 | (g) 



, (5) 



dt 15 M2 p4 \^ ' 
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In these equations /i = 17111712/ {mi +777.2) is the reduced mass and, mi and m2 
are the two component masses. We can convert Eqs. (4), (5), and (6) into the same 
notation as in [3], by replacing p by a using Eq. (2). 

In order to compute eccentric waveforms, we need to solve the system of 3 ODEs 
as defined by Eqs (4), (5), and (6). The method employed is described in details in 
Sec. 2.2. 

Once the evolution of the three parameters e, p and (j) are known, we can generate 
the eccentric waveforms. The source and the gravitational wave detector are separated 
by a distance, R. The detector is placed in a direction defined by the polar angles i 
and P relative to the Cartesian frame (following standard convention [16]). The unit 
vectors defined by t and (3 are chosen as polarization axes. The following equations 
are the two fundamental polarizations of the gravitational waves [13]: 



<5+ 



5e 

2 cos(2(/) - 2/3) + y cos((/) - 2/3) 



■ cos(30 - 2/3) + cos(2/3) 



(1 + cos^ l) 



[e cos(0) ■ 



e^] sin^ L 



(7) 



-4^-4 sin(2(^ - 2/3) + be sin((^ - 2/3) 
pR I 

esin(3(^-2/3) - 2e^ sin(2/3) cos i. 



(8) 



In Sec. 3.2, we will see how t and /3 affect the waveform shape. It is important to 
notice that the signal can be decomposed into components that oscillate at once, twice 
and three times the orbital frequency. This effect is entirely related to the presence of 
eccentricity. Indeed if we set cq = 0, the remaining terms in Eqs (7) and (8) contain 
components that oscillate at twice the orbital frequency only. 

We also need the initial conditions for the semi-latus rectum, po, the orbital phase, 
(/)o, and the eccentricity, cq. The initial conditions on cq and are set arbitrarily and 
are specified in Sec. 3. 2. The remaining initial condition, po^ is derived from the initial 
eccentricity, cq. Using Eq. (3), we have 



1 



Po 



(9) 



{27rMfof' ' 

where /o depends on the lower cut-off frequency, F^^ at which the signal enters a 
detector's sensitivity band. If we want the waveform to be valid from Fl onwards, 
we must have the third harmonic to start at Fl. The consequence is that the first 
harmonic, which contributes to the third, must start at Fl/3. So, in Eq. (9), we must 
set /o = Fl/3. 

Finally, let us have a short digression about the effect of the parameter /o, which 
may lead to difficulties from the point of view of data analysis. Indeed, if eo = then 
the first and third harmonic contributions are null (see Eqs. (7) and (8)). Yet, with the 
definitions provided above, we can see that the second harmonic starts at 2Fl/3. In 
the case of a non-eccentric waveform implementation, starting the waveform at Fl is 
sufficient. So, if we compare the duration of an eccentric waveform (even when eo = 0) 
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as defined in ttiis paper witii tiie duration of a purely non-eccentric waveform then the 
ratio of the durations of two waveforms is (2/3) For instance, a (1.4 + 1A)Mq 
binary lasts 73 seconds if Fl = 40 Hz, and 465 seconds if Fl = 20 Hz. These 
durations have to be compared to a standard implementation where = Fl : 25 and 
158 seconds, respectively. Dealing with such long waveforms may lead to technical 
issues from a computational point of view (e.g., memory allocation). 

2.2. Eccentric waveforms: validation and investigations 

The main difficulty in generating eccentric waveforms from Eqs. (7) and (8) resides 
in the resolution of the system of differential equations. This is done by using a 
numerical code based on the GNU Scientific Library [19]. The ODE integrator is 
based on a Runge-Kutta-Fehlberg (4, 5) method. The step size is taken to be the 
inverse of the sampling frequency (4 kHz). Using this method, we compute the 
waveform durations as a function of initial eccentricity and component masses (see 
Table. 1). Our results do not agree with those obtained in [11]: Durations generated 
by our code are consistently 5% lower than those computed in [11]. This effect has 
been seen independently [27]. We performed another sanity check in the case where 
eo = 0: we compare our waveform with a non-eccentric waveform (at Newtonian 
order) that is generated with the LIGO Algorithm Library (LAL) [18]. We found a 
perfect agreement^ between the two waveforms. So, we think that the results in [11] 
may be biased due to a lack of accuracy in the evolution of the ODEs. We are now 
able to extend the waveform generation to higher masses and larger eccentricities as 
compared to Ref. [11], whose study was limited to (8.0 + 8.0) M©. 

The discrepancy that has been noted and the ability to extend investigations to 
higher mass range is also one of the motivation for revisiting their work. 

In Table. 1, we provide the duration of the eccentric waveforms, denoted 
T^cq^Fl)^ for various component masses and eccentricities. Table 1 shows that 
for a given total mass, the duration of the signal decreases with increasing initial 
eccentricity. If a binary is highly eccentric initially, then its orbit shrinks faster by 
losing energy at a faster rate via radiation of gravitational waves. Moreover, for a 
given initial eccentricity, the waveform duration decreases as the total mass increases. 
This behavior can be anticipated by looking at Eq. (5), for which a larger initial value 
is expected for larger masses and larger eccentricity. 

Using Tab. 1, we can also compare the duration of the circularized waveform, 
^c(eo = 0), to the duration of the eccentric waveform, T(eo). The ratio of these 
two durations is shown in brackets. It is interesting to notice that the ratio 
T(eo)/Tc(eo = 0) depends on the eccentricity only and is independent of the total 
mass. Knowing the duration of the circularized waveform, we can anticipate the 
duration of the eccentric waveform by using the following expression 

T(eo) = Tc(eo = 0) [l + 0.02eo - 4.37eg + 3.53eg] . (10) 

Although the derivation of Eq. (10) is empirical, it is accurate to 1% error for 
eccentricity up to 0.6. This expression is independent of the initial cut-off frequency 
and masses of the system (in the Newtonian case) as shown in Fig. 1. This expression 
could be useful for data-analysis purposes so as to optimize the vector's length to be 
stored in memory. 

X To obtain the same waveform duration, we need to set the lower cut-off frequency of the non- 
eccentric waveform to 2Fl/?> Hz. 
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eo 


1.4+1.4 


1.4+10 


5+5 


10+10 


20+20 





73.30(1.00) 


16.37(1.00) 


8.77(1.00) 


2.75(1.00) 


0.85(1.00) 


0.1 


70.66(0.96) 


15.78(0.96) 


8.46(0.96) 


2.65(0.96) 


0.82(0.96) 


0.2 


63.18(0.86) 


14.11(0.86) 


7.56(0.86) 


2.37(0.86) 


0.73(0.86) 


0.3 


52.09(0.71) 


11.63(0.71) 


6.23(0.71) 


1.95(0.71) 


0.60(0.71) 


0.4 


39.13(0.53) 


8.73(0.53) 


4.68(0.53) 


1.46(0.53) 


0.43(0.50) 


0.5 


26.22(0.36) 


5.84(0.36) 


3.13(0.36) 


0.97(0.35) 


0.27(0.31) 


0.6 


15.09(0.21) 


3.35(0.20) 


1.80(0.20) 


0.55(0.20) 


0.06(0.07) 


0.7 


6.92(0.09) 


1.52(0.09) 


0.82(0.09) 


0.23(0.08) 


/ 


0.8 


2.12(0.03) 


0.40(0.02) 


0.23(0.02) 


/ 


/ 



Table 1. Duration of the gravitational- wave signal as a function of initial 
eccentricity and component masses. The lower cut-off frequency is Fl = 40 Hz. 
The numbers provided in this table show a systematic difference of about 5% 
compared to numbers from Tab. I of [11]; we believe that our results are 
more accurate than those from [11] (see the text for an explanation). Because 
the duration decreases when eccentricity and total mass increase, there is a 
combination for which the ending frequency of the waveform is below F^. In 
such a case, no waveform can be generated, which is represented by the / symbol. 
In brackets, we also provide the ratio of the duration of an eccentric waveform 
with respect to the first line, where eo = 0. This ratio is independent of the 
component masses (See Fig. 1 and the text for more details). 
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Figure 1. The ratio of the duration of an eccentric binary to that of a quasi- 
circularized waveform as a function of eccentricity. The durations of the eccentric 
(T) and circularized (Tc) systems are extracted from the integration of the ODEs. 
The ratio as a function of eccentricity can be accurately fitted with a cubic 
equation, which depends only on the initial eccentricity cq. The crosses and 
circles show the ratio T/Tc for two different lower cut-off frequencies F^ = 20 Hz 
and Fl = 40 Hz. The fit is also independent of F^. It is also independent of the 
total mass (See Tab. 1). 

Let us finish this section with an example. Let us consider a (5.0 + 5.0) Mq 
system with cq = 0.5. Duration of such a waveform is 3.13 s. We want to emphasize 
that the parameters /3 and i are of no importance for the simulations we present in 
Sec. 3. In Fig. 2, we plot the first 0.5 s of this waveform when P = and l = ir/A. 
Then, we keep (3 = and set ^ = 0. We see that only the amplitude of the waveform 
changes; the period remaining unchanged. Finally, we set /3 = 7r/4 and l = tt/A. Now, 
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the amplitude changes as weh as the phase. However, the phase is only shifted. These 
results are expected if we look at Eq. (7) and Eq. (8), where we see that /3 contributes 
to the phase, (/), in the same manner in all the three harmonics, and t comes as a factor 
for the amplitude only. So, neither f3 nor t affect the overall frequency behavior of the 
waveform. 




Figure 2. Variation of a (5.0 + 5.0) Mq eccentric waveform with an initial 
eccentricity eo =0.5 as a function of different (3 and l values. 



3. Searching for eccentric binaries with templates from non-eccentric 
binaries 

In this section, we briefly describe the different tools that we have used in our 
simulations. Then, we present the simulation protocol that we have followed. Finally, 
we present the results of our simulations in the context of different detectors and mass 
ranges. 



3.1. Match and template bank 

We denote the Fourier transform of two functions x{t) and h{t) by x{f) and h{f), 
respectively. The matched-flltering inner product is deflned by 



i 



°° x*{m f)+xi f)h*if) 

Skif) 



ix,h) = 4 ^MJl^MTlJ^yUllMldf , (11) 



where Sh{f) is the noise power spectral density (PSD) of the detector. A normalized 
signal is given by 

^^^TP^' ^^^^ 

where a hat denotes normalized signal. The SNR is deflned by 

''"=Jn5 = (''')^ '''' 



We can now imagine that x{t) = n{t) + s{t^i9^), where n{t) is Gaussian noise 
and s{t^i^^) is a gravitational-wave signal characterized by a set of p parameters 
'i^^ , /i = 0, 1, — 1. These parameters may be separated into intrinsic and extrinsic 
parameters. The intrinsic parameters are represented by the component masses mi 
and m2 and the initial eccentricity eo, while the extrinsic parameters are the initial 
orbital phase 0o and the time of arrival, to. When the noise fluctuations are neglected, 
we can define the match between a signal s{t, "d^) and a template h{t^ as follows [25]: 

M {s{t, d^), h{t, ^^)) = max (s{t, d^), h{t, ^^)) (14) 

^0,00 V / 

where the extrinsic parameters are maximized over automatically (using Fourier 
transform). 

Using a continuous template space, we can maximize the matches over the 
intrinsic parameters to obtain the fitting factor, FF, between the signal and the 
template family [28], which is defined by 

FF {s(t, d^) , h(t, r)) = max M (s(t, d^), h(t, i?^)) . (15) 

The fitting factor can be interpreted as the maximum fraction of SNR that can be 
obtained by filtering a signal with an approximate template family. The signal and 
template families can be different, like in this paper, where the signal parameters are 
'i^^ = mi, m2, eo and the template parameters are i^^ = mi^m2. The fitting factor was 
used in Ref. [11] to estimate the loss of SNR when an eccentric signal is filtered with 
a non-eccentric signal. 

In practice, the continuous template space is replaced by a discrete one, which 
is called a template bank [25, 22, 23]. It is represented by {h{t^i9^)}^ where 
i = 0,1,..,A^5 — 1, and is the number of templates. A template bank is optimally 
designed if N^j is the smallest such that for any signal there always exists at least one 
template in the bank that gives 

minmaxM(5(t,i?^),/i(t,i?n) > MM, (16) 

where MM is the minimal match defined by the user. Usually, MM is set to 95%, or 
97%, which corresponds to a loss in detection rate of about 14% and 9%, respectively 
(the expected loss in detection equals 1 — MM^). So, the minimal match is a measure 
of how well a discrete template bank covers the parameter space defined by the intrinsic 
parameters. 

Let us define a new notation that will be useful to quantify our results. We define 
the match over the bank, M^, as follows 

MB{s{t, d^), h{t, = max M (^s{t, d^), h(t, i^^) • (17) 

The quantity 1 — Mb can be interpreted as the average fractional loss of SNR. 

Searches for inspiralling compact binaries performed in LIGO data use a template 
bank constructed so that signal from any circularized binary, s{t), is found with a 
match greater than the minimal match, MM. Exhaustive simulations were performed 
in [26, 20] to test the template bank placement in various cases: mass parameters 
related to BNS, BBH and BHNS systems and various design sensitivity curves of 
ground-based detectors. 

In the rest of this paper, we will study the distribution of the quantity Mb when 
eccentric binaries are searched with the same template bank placement as the one 
used in LIGO analysis, where templates are generated with a non eccentric model. 
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initial LIGO 



Adv. LIGO 



E.T. 



VIRGO 



Bank size 
Bank size BNS 



6792 
2319 



25194 
6969 



65104 
16046 



52455 
12955 



Table 2. Number of templates, A^^, of each template bank used in our simulations. 
The template bank size of initial LIGO is the smallest because its lower cut-off 
frequency, Fl equals 40 Hz whereas other detectors have Fl = 20 Hz. The 
parameter space of the first row is mi,m2 G [1,3O]M0 for initial LIGO and 
TTii^ 1712 G [1, 4O]M0 for the others. The second row gives A^^ for a BNS parameter 
space, where mi,m2 G [1,3]M0 (see the appendix). 



3.2. Simulation parameters 

As mentioned in the introduction, we want to identify the parameter space (in total 
mass, M, and initial eccentricity, eo) where searching for eccentric binaries with 
templates corresponding to circularized systems suffices to obtain a negligible loss 



Our signals are based on eccentric waveforms whose phase is at Newtonian order 
only, we will therefore also consider Newtonian order for our templates. We will use 
the so-called TaylorTS approximant (see e.g, [20, 18] for a precise definition). The 
lower cut-off frequency of the template is chosen to be Fl. 

In all the simulations that follows, we use an hexagonal template bank with 
a minimal match MM = 95%. It has been shown that for VIRGO, LIGO, Einstein 
Telescope and advanced LIGO design sensitivity curves, Mb is guaranteed to be larger 
than the minimal match (95%) for any system with 2Mq < M < 60Mq [20]. Yet, 
in some cases, as described later in Sec. 3.3, we will extend the search to 8OM0. 
Note that the template bank is actually optimized for 2PN order in phase and might 
therefore not be optimal for the Newtonian order considered here. In principle, since 
the phase of the signal and template are based on a Newtonian order only, we could use 
a one dimensional template bank (e.g., the chirp mass as defined later). So, using the 
template bank described in [20] and used in LIGO searches, which is a two dimensional 
grid, our simulations may provide slightly over estimated values of matches. 

In all our simulations, we use the following common parameters. The number 
of waveforms in each simulation is 10,000. The sampling frequency is 4096 Hz. The 
total mass is uniformly distributed in the range considered (see Sec. 3.3). The initial 
starting phase is randomized between and 27r. Although the choice of the initial 
phase should be of no consequence, we randomize it in all our simulations. The initial 
eccentricity is uniformly distributed in the range considered, that is eo G [0.0, 0.4]. We 
also fix^ = 7r/4, /3 = 0. As explained in Sec. 2.2, fixing P and i to arbitrary values 
does not affect the waveform significantly. However, we performed tests where p and 
L were randomized. As expected, we did not see any significant effects on the matches. 
Finally, let us note that if f3 varies, then the ending frequency may slightly change 
and the matches as well. 

3.3. Search for eccentric waveforms in ground-based detectors 

We consider 4 design sensitivity curves: initial LIGO, advanced LIGO, Einstein 
Telescope and VIRGO. In this section we summarize the results we obtained when 
searching for eccentric waveforms with circular waveforms. We limited the maximum 
value of eccentricity to eo = 0.4, above which the loss of SNR is always larger than 
50%. Note that for all the results presented below, the initial eccentricity is defined by 



of SNR. 
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Eq. (3), which is defined at /min; in other words, when the main harmonic frequency 
equals ^Fj^jZ. 

S.S.I. Initial LIGO In the case of initial LIGO, the lower cut-off frequency is 
Fl = 40 Hz and the total mass range M G [2,6O]M0. The template bank size is 
provided in Tab. 2. In Fig. 3, we show the distribution of the match over the bank, 
Mb, in the (M, eo) plane. As expected, around eo = 0, Mb is close to 1. However, 
in most of the parameter space, Mb is much smaller than unity. It is convenient to 
represent Mb with isocontours set to [96.5, 95, 90, 80, 50]%. Using a value of 96.5% is 
useful because the loss of detection associated with it is about 10%. Using a value of 
95% is also useful because it corresponds to the minimal match of the template bank. 




Eccentricity 

Figure 3. Maximum matches over the template bank (Mb) between an eccentric 
signal defined by its total mass (y-axis) and initial eccentricity at 2Fl/3 = 
26.66 Hz (x-axis). The design sensitivity curve is based on initial LIGO. The colors 
represent the value of the maximum match over the template bank. The minimal 
match of the bank is 95%. So the area where matches are below 95% represent 
the region of parameter space where eccentricity significantly affect search for 
eccentric binaries with template based on circularized system. 

The isocontours follow more or less the same structure. If we focus on the 
95% isocontour, we can see that BNSs are found with matches above 95% if initial 
eccentricity is below 0.05. While the total mass increases up to 20 M©, the initial 
eccentricity for which the match is greater than 95% increases up to a maximum value 
of 0.15. As the total mass increases for a match of 95% or greater, from 20 Mq to 
about 40 M©, the initial eccentricity again goes down to 0.05. Above 40 M©, matches 
can go down to 90% even for negligible eccentricities but this is related to the template 
bank design. Indeed, the template bank is optimally designed for signals that have a 
significant part of their power in a frequency band that resides where the detector has 
the best sensitivity. This assumption breaks down for signals with M > 4OM0, which 
ending frequency is already as low as 110 Hz whereas the detector best sensitivity is 
about 200 Hz. Consequently the bank is sub-optimal for these systems. 

In the appendix, in Fig. Al, we provide for convenience the same results as 
in Fig. 3 where M has been replaced by the chirp mass, A^, which is defined by 
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M = Mrf^^ and r] = ^^^^ . We also provide in Fig. A2 the same results as in Fig. 3 
where Mb has been replaced by the loss of detection (1 — M^). 

Finally, since the parameter space related to BNSs is rather small in Fig 3, we 
performed another simulation where we focus on the range M G [2,6]M0. As we can 
see in Fig A3, using non-eccentric waveforms is sufficient to detect a (1.4 + 1A)Mq 
eccentric waveforms if eo < 0.05. 

3.3.2. Advanced LI GO The simulation in the case of advanced LIGO differs from 
that of initial LIGO in two ways: the lower cutoff frequency is set to Fl = 20 Hz, and 
the maximum total mass is extended to 80 Mq. Results are shown in Fig. 4 and are 
very similar to the initial LIGO case. The main difference being that the high mass 
range is extended because the lower cut-off frequency is smaller. 

If we focus on the 95% isocontour, we can see that BNSs are found with matches 
above 95% if initial eccentricity is below 0.05. Then, for the same match the initial 
eccentricity increases up to a maximum value of 0.15 as the total mass increases up 
to 25 Mq. Finally, when the total mass increases from 25 Mq to about 70 M©, the 
initial eccentricity remains constant around 0.15 for the same match of 95%. Above 
70 M0, matches go down to 90% even for negligible eccentricities. 

In the appendix (Fig. Al), we provide for convenience the same results as in Fig. 3 
where M has been replaced by the chirp mass. In Fig. A2), we also provide the same 
results as in Fig. 3 where Mb has been replaced by the loss of detection (1 — M^). We 
performed another simulation so as to focus on the ENS region. Results are shown in 
Fig. A3: using the non-eccentric model is sufficient to detect (1.4 + 1.4)M0 eccentric 
systems with eo < 0.05. 




0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 

Eccentricity 



Figure 4. Same as Fig. 3 but for advanced LIGO. 



3.3.3. Einstein Telescope The simulation in the case of the Einstein Telescope is 
similar to that of advanced LIGO (same F^^ same mass range). Results are shown in 
Fig. 5. 

If we focus on the 95% isocontour, we can see that BNSs are found with matches 
above 95% if e ^ 0. Then, the initial eccentricity increases up to a maximum value of 
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0.10 while total mass increases up to 30 Mq while obtaining a match ^95%. Finally, 
when the total mass increases from 30 Mq to about 60 M0, the initial eccentricity 
decreases to eo ~ 0. Among the four design sensitivity curves investigated in this 
paper, the Einstein Telescope case gives the least satisfactory results. 

In the appendix, in Fig. Al, we replace M by using results from Fig. 5. In 
Fig. A2, we replace Mb by the loss of detection (1 — M^). Finally, we focus on the area 
M e [2,6]M0 and provide the results in Fig A3: using non-eccentric model sufficient 
to detect a (1.4 + 1.4)M0 eccentric system with cq < 0.02. 




Eccentricity 

Figure 5. Same as Fig. 4 but for the Einstein Telescope. 



3.3.4' VIRGO The simulation in the case of VIRGO is similar to that of advanced 
LIGO (same Fl^ same mass range). The results are shown in Fig. 6. 

If we focus on the 95% isocontour, we can see that BNSs are found with matches 
above 95% if initial eccentricity is below 0.03. Then, the initial eccentricity can be 
increased up to maximum value of 0.15 while total mass increases up to 40 Mq which 
still achieve matches ^ 95%. Finally, when the total mass increases from 40 Mq to 
about 80 M0, the initial eccentricity decreases to cq = 0.15 for matches ^ 95%. 

In the appendix, in Fig. Al, we replace M by using results from Fig. 6. In 
Fig. A2, we replace Mb by the loss of detection (1 — M^). Finally, we focus on the 
area M G [2^6]Mq and provide the results in Fig A3: using non-eccentric model is 
sufficient to detect a (1.4 + 1.4)M0 eccentric binaries with cq < 0.03. 

4. Conclusion and perspectives 

In this paper the main difference with previous works is that we maximize the matches 
over a discrete template bank, M^, instead of computing the so-called fitting factor, 
FF, which maximizes the matches over a continuous parameter space. The template 
bank we used is identical to the template bank used in LIGO searches for binaries in 
quasi-circular orbits. For this reason, and since our estimation is more accurate, our 
results are slightly different from [11]. 
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Figure 6. Same as Fig. 4 but for VIRGO. 

If we search for gravitational-wave signal emitted by eccentric inspiralling 
compact binaries with standard non-eccentric template bank, our main results can 
be summarized as follows. In initial LIGO, (1) any binary neutron star system can 
be detected with Mb > 95% if the initial eccentricity eo < 0.05 (or Mb > 90% if 
eo ^ 0.10) and (2) any binary system with a total mass, 6 < M < 35 Mq can be 
detected with Mb ^ 95% if eo ^ 0.10. Similar results have been obtained for the 
advanced LIGO case. Concerning VIRGO detector, any binary neutron star system 
can be detected with Mb > 95% if eo < 0.04, or Mb > 90% if eo < 0.06, and 
(2) any system with total mass larger than 6 Mq can be detected with Mb ^ 95% 
if ^0 ^ 0.10. Similar results have been obtained in the case of Einstein Telescope. 
In conclusion, the initial eccentricity does not need to be taken into account in our 
searches for gravitational- wave signals if eo ^ 0.05, because the loss of SNR (about 5%) 
is consistent with the discreteness of the template bank that is being used. However, 
by neglecting the eccentricity in our searches, significant loss of detection may arise as 
soon as eo > 0.05, especially for binary neutron stars and systems with high masses 
(M > 35M0). 

Future directions for this work are the inclusion of higher PN order, the 
development of a template bank placement for eccentric systems and parameter 
estimation. Concerning higher PN order, we need to reproduce the simulations 
provided in this work. This work should be straightforward given the infrastructure 
which has been developed in this paper. Concerning the template bank placement, 
as we have explained, above eo = 0.05 the loss of SNR may be quite large in some 
region of the parameter space. A solution may be to use a dedicated template bank 
to take into account the eccentricity. For instance, using the template bank used in 
this paper, we can have several layers of templates, each having a different values of 
initial eccentricity. More investigations are needed but preliminary studies shows that 
only a few layers are required to significantly improve the results (the TaylorT3 family 
being replaced by the eccentric waveform presented here). Finally, using Monte-Carlo 
simulations it would be interesting to estimate the error bound on the variance of the 
measured eccentricity. 
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Appendix A. Additional results 

In Fig. Al we provide the results shown in Sec. 3. 2 in term of chirp mass rather than 
the total mass. See the text for a full description of the simulation parameters. We 
also transform the results into loss of detection in Fig. A2. Finally, we further explore 
BNS region by carrying out a simulation restricted to the mass range M G [2,6]M0, 
and the results of which are shown in Fig. A3. 
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Figure Al. Matches between the simulated eccentric waveform and a template 
bank made of non-eccentric waveforms. The eccentricity is uniformly distributed 
in the range [0,0.4]. The y-axis represent the chirpmass parameter. The data used 
are identical to those of Fig. 3, 4, 5 and 6, where the design sensitivity curves are 
from top left to bottom right : initial LIGO, advanced LIGO, Einstein Telescope 
and VIRGO. 
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Figure A2. Loss of detection. Simulated waveform have eccentricity in the 
range [0, 0.4]. Template bank is made of non-eccentric waveforms. From top left 
to bottom right, we used initial LIGO, advanced LIGO, Einstein Telescope and 
VIRGO design sensitivity curves. The data used are identical to those of Fig. Al 
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Figure A3. Matches between the simulated eccentric waveform (BNS) and 
a template bank made of non-eccentric waveforms. The initial eccentricity is 
uniformly distributed in the range [0, 0.1]. From top left to bottom right, we used 
initial LIGO, advanced LIGO, Einstein Telescope and VIRGO design sensitivity 
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